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^ ' We determine the second order endpoint of the line of first order phase transitions, 
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and analyze universal properties of this chiral critical point. A detailed analysis of 
Binder cumulants and the joint probability distributions of energy like and ordering- 
field like observables confirms that the chiral critical point belongs to the universal- 
ity class of the three dimensional Ising model. From a calculation with improved 
^' gauge and staggered fermion actions we estimate that the transition is first order 
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for pseudo-scalar meson masses less than about 200 MeV. 
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1 Introduction 



In the chiral limit the order of the QCD phase transition depends on the number 
of quark flavours that become massless. It has been conjectured [jl| and verified 
in numerical calculations 0, |^, ^ that this phase transition is first order for QCD 
with Uf > 3 massless quark flavours. The strength of this first order transition 
weakens in the presence of non- vanishing quark masses and eventually ends in a 
second order endpoint. The analysis of effective models constructed in terms of the 
chiral order parameter suggests that this endpoint, the chiral critical point, belongs 
to the universality class of the three dimensional Ising model |^. In fact, the entire 
surface of second order phase transitions which in the n/-dimensional space of quark 
masses separates the regime of first order transitions at small quark masses from 
the crossover region is expected to belong to this universality class. 

In the case of degenerate quark masses the lattice formulation of QCD depends 
on two bare couplings, the bare quark mass m and the gauge coupling (3 = 6/g^- The 
bare quark mass acts like an external symmetry breaking field and thus leads to an 
explicit breaking of the SUL{nf) x SUR{nf) chiral symmetry. For nj > 3 there is a 
first order phase transition in the chiral limit (m = 0), which will continue to persist 
for small but non- zero values of m up to a critical value, m, of the quark mass. It is 
expected that the universal properties of this chiral critical point are controlled by 
a global Z{2) symmetry, which however is not an obvious global symmetry of the 
QCD Lagrangian. It rather is the relevant symmetry of the effective Hamiltonian 
which controls the critical behaviour at this point. Although the transition is second 
order at the chiral critical point, it is obvious that neither the chiral condensate nor 
the Polyakov loop will be an order parameter for the spontaneous Z(2) symmetry 
breaking at this critical point. In fact, we do not know a priori what are the 
relevant observables representing the energy-like and ordering-field like operators of 
the effective Hamiltonian and the corresponding scaling fields (couplings) that define 
the energy-like and ordering-field like directions in the coupling parameter space; it 
is part of the problem of analyzing the universal behaviour at this endpoint that 
one has to identify these operators and couplings appropriately. 

The problem of determining the critical properties at the endpoint of a line of 
first order phase transitions is well known from studies of other statistical and field 
theoretic models. The construction of appropriate scaling fields has been discussed 
in detail in the case of the liquid-gas phase transition 0. The concepts developed 
in this context have recently also been used to locate and explore the properties of 
the critical endpoint of the electro- weak phase transition |^, Q as well as the critical 
point in the ferromagnetic, three dimensional 3-state Potts model ^. The latter 
problem, of course, is closely related to the analysis of the endpoint of the line of 
first order phase transitions that occur in the heavy quark mass limit of QCD |1C]. 



We will focus here on a discussion of the light quark mass regime of QCD. 
In particular we concentrate on an analysis of universal properties of the chiral 
critical point in 3-flavour QCD. Although we will eventually also be interested in 
an accurate determination of the location of this endpoint in terms of physical mass 
values, e.g. quark or hadron masses, our main emphasis here will be on the analysis 
of universal properties. We thus will mainly use standard (unimproved) Wilson 
gauge and staggered fermion actions for our analysis. Lattice calculations will be 
performed on lattices with temporal extent N^ = 4 and various spatial lattice sizes 
to use well known finite size scaling methods at the critical point. A first estimate 
for the location of the critical endpoint on lattices with temporal extent A^^ = 4 has 
been given in |Tl|. We confirm this estimate here and will, in addition, present a 
detailed analysis of the universal properties of the transition at the chiral critical 
point. Using improved gauge and staggered fermion actions we then will give a first 
estimate of the pseudo-scalar meson mass at the chiral critical point. 

This paper is organized as follows. In the next section we introduce the basic 
variables relevant for the description of critical behaviour in the vicinity of the 
chiral critical point. In section 3 we locate the chiral critical point and determine 
its universality class using Binder cumulants. In section 4 we construct the order 
parameter and energy-like operator at the chiral critical point. Section 5 is devoted 
to the determination of the critical pseudo-scalar meson mass at the critical point. 
Finally we give our conclusions in Section 6. 



2 Scaling fields at the chiral critical point 



The thermodynamics of QCD with Uf quark flavours is described in terms of the 
partition function 



Z(/5,m)= / pf/e-^"^}'/^''") 



which is defined on a 4-dimensional lattice of size A^^ x A^. Here (3 = Q/g"^ and 
m denote the gauge coupling and bare quark mass, respectively. The Euclidean 
action, S, is given in terms of the pure gauge action Sg and the fermion matrix Qp- 
Throughout this paper we will use staggered fermion actions and standard Wilson 
or tree level improved gauge actions. The discretized QCD action then reads. 



S{{U},P,m)=pSG{{U})-'^TT lng^({f/},m). (2) 




Figure 1: Generic phase diagram in the low mass regime of 3-flavour QCD. 

The generic phase diagram of QCD in the small quark mass region is shown in 
Fig. |I| in the plane of bare couplings appearing in the QCD partition function. In the 
vicinity of the second order endpoint, {(3c{fn) , fn) , the dynamics and the universal 
critical behaviour is controlled by an effective Hamiltonian, which can be expressed 
in terms of two operators S, Ai, i.e. the energy-like and ordering-field like operators 
that couple to two relevant scaling fields r and ^, 



He//(r,0 = r£: + e-M 



(3) 



Under renormalization group transformations the multiplicative rescaling of the cou- 
plings ^ and T is controlled by the two relevant eigenvalues that characterize the 
universal critical behaviour in the vicinity of the second order critical point. The 
singular part of the free energy density thus scales like 



fs{T,o = b-'fs{byw,by^O 



(4) 



where the dimensionless scale factor b = LT = N„/Nj. gives the spatial extent of 
the lattices in units of the inverse temperature. Susceptibilities constructed from S 
and Ai will show the standard finite size scaling behaviour 



X£^V-\{6£f)r^asb^/'' , XM^V-'{{6Mf)^aMb''/'^ 



(5) 



Here V = N^ denotes the 3- dimensional spatial volume and 6X = X — (X). Note 
that S and A4 as well as Sq and ipip are defined to be extensive quantities. 

For finite values of the quark mass all global symmetries of the QCD Lagrangian 
are explicitly broken. As the symmetry of Tie// that characterizes the critical be- 
haviour at the chiral critical point is not shared in any obvious way by the QCD 
Lagrangian we also may expect that in the vicinity of the chiral critical point the 
operators appearing in the QCD Lagrangian are mixtures of the energy-like {S) and 
ordering-field like [Ai) operators. Similarly the couplings appearing in the QCD 
Lagrangian are linear combinations of the scaling fields as we have indicated in 

Fig. m. 

In the vicinity of the critical point one may use a linear ansatz for the couplings 



C = m-m + B{p-p,) , (6) 



as well as for S and Ai which are constructed in terms of operators appearing in 
the original QCD Lagrangian, 



M = ipij + sSc . (7) 



Here ipip denotes the chiral condensate evaluated on a given gauge field configuration 
{U}. In terms of the staggered fermion matrix, Qp{m), this is given by ipip = 
0.25 nfTTQp\{U},m). 

As the operators of the QCD Lagrangian, e.g. Sg and ipip or related observables 
like the Polyakov loop expectation value, are mixtures of S and Ai, the correspond- 
ing susceptibilities will all receive contributions from fiuctuations of S as well as 
A4. Asymptotically therefore all of them will show identical finite size scaling be- 
haviour which will be dominated by the larger of the two exponents a/u and 7/1^, 
respectively. For the symmetry groups of interest in the QCD context, e.g. the 
symmetry of three dimensional Z{2) or 0{N) spin models, this will be 7/z/. A finite 
size scaling analysis of susceptibilities constructed from the basic operators of the 
QCD Lagrangian thus will give access only to the ratio 7/?^, which unfortunately is 
quite similar for all the above mentioned symmetry groups and thus is not a good 
indicator for the universality class controlling the critical behaviour in the vicinity 
of the chiral critical point. 



The situation is different for cumulants constructed from linear combinations of 

ijjip and Sc, 



{(6M(x)Y) , , - 



((<5M(x))^) 



2 



From Eq. ^ it follows that for arbitrary values of x the cumulants are renormalization 
group invariants which in the infinite volume limit take on a universal value at the 
critical point (r, ^) = (0,0). For all values of x different from 1/r the cumulants 
behave asymptotically like the Binder cumulant for the order parameter; cumulants 
calculated on different size lattices for different quark masses will intersect at some 
value of the quark mass. In the infinite volume limit these intersection points will 
converge to a universal value which is characteristic for the universality class of the 
underlying effective Hamiltonian and, in fact, is quite different for the classes of 
three dimensional 2'(2) and 0{N) symmetric spin models; e.g. S4 = 1.604 for Z{2) 
g, 1.242(2) for 0(2) |13[ and 1.092(3) for 0(4) H. The cumulants B^^x) thus 



seem to be appropriate observables to locate the chiral critical point as well as to 
determine its universality class without knowing in detail the correct scaling fields. 



3 Locating the chiral critical point 



The determination of the chiral critical point proceeds in two steps. First of all, 
we determine for fixed values of the quark mass pseudo-critical couplings, Ppc{m), 
on finite lattices. These are defined as the position of maxima in susceptibilities of 
ipip, the gauge action Sq and the Polyakov loop L [jl5|. We then make use of the 
finite size scaling properties of Binder cumulants -84(0;) evaluated at (3pc{m). When 
analyzed as function of the bare quark masses the cumulants calculated on lattices 
with spatial extent Li and L2 will intersect at a mass mLi,L2- For {Li, L2) -^ (00, 00) 
these intersection points will converge to the chiral critical point. 

From previous studies with standard Wilson gauge and staggered fermion actions 
one knows that the endpoint in 3-flavour QCD is located close to m = 0.035 |Tl| . 



As the universal properties of the endpoint are not expected to be influenced by 
lattice cut-off effects we took advantage of this knowledge and performed our detailed 
scaling analysis with unimproved actions. We have performed calculations on lattices 
of size Nl X 4, with N^j = 8, 12 and 16. We have used four values of the quark 
mass in the interval m G [0.03, 0.04] and for each of these masses we calculated 
thermodynamic observables for 3 to 4 different values of the gauge coupling (3. In 
general we collected for each pair of couplings (1 — 3) ■ 10^ configurations generated 

6 



with the hybrid-R algorithm^. Interpolations between results from different /3-values 
have been performed using the Ferrenberg-Swendsen multi-histogram technique |TB . 



m 


iV. 


v-^ {{5^i^f) 


V-' {{6Lf) 


V-' {{SSgT) 


0.0300 


8 


13.4(2) 


1.32(3) 


3.21(10) 




12 


32.5(1.4) 


3.08(12) 


7.34(31) 




16 


65.0(4.7) 


6.27(44) 


14.8(1.0) 


0.0325 


8 


12.9(4) 


1.32(4) 


3.21(12) 




16 


50.5(3.0) 


5.00(29) 


11.72(70) 


0.0350 


8 


11.9(5) 


1.34(5) 


3.05(14) 




12 


23.5(1.5) 


2.52(16) 


5.73(40) 




16 


39.8(2.8) 


4.19(28) 


9.52(64) 


0.0400 


12 


16.8(8) 


2.05(11) 


4.45(24) 




16 


23.4(2.1) 


2.84(23) 


6.03(54) 



Table 1: Volume dependence of the susceptibilities 



In Tab. |I| we summarize our results for the peak heights of the three different 
susceptibilities which we have analyzed. We generally find that the positions of 
these peaks coincide within statistical errors. Thus in Tab. |^ only the pseudo- 
critical couplings extracted from the location of the peak in {{dipipY) are given. 
Also given in this table is the Binder cumulant B/i{x) calculated at the pseudo- 
critical couplings for two values of x. We have checked that the cumulants indeed 
attain their minimum at (3, 



pc 



rn). The case x = corresponds to the cumulant of 
the chiral condensate alone, whereas x = 0.43 corresponds to our best estimate for 
the mixing parameter s, whose determination we are going to discuss in the next 
section. 

The cumulant of the chiral condensate, -04(0), is shown in Fig. ^. We note that 
the cumulants calculated on different size lattices intersect at a quark mass close 
to m = 0.035. The value of -64(0) at the intersection point is compatible with the 
universal value of the Binder cumulant for the 3-dimensional Ising model. It is 
obvious from the data given in Tab. ^ that the situation is very similar for x = 0.43. 
We have fitted the cumulants on a given lattice size using a linear ansatz in the 
quark mass, B^lx) = ao + cti m. From this we find for the intersection point and 
the cumulant, 



X 



X 



= 0.0 
0.43 






m] 



m] 



^We used trajectories of length r 



(5.1458(5), 0.0331(12)) , ^4(0) = 1.639 ± 0.024 , 
(5.1454(5), 0.0329(15)) , B^{x) = 1.624 ± 0.023 . (9) 

0.675 generated with a discrete step size 6t — 0.015. 



m 


N. 


Ppc 


^4(0) 54(0.43) 


0.0300 


8 

12 

16 


5.1403(5) 
5.1411(4) 
5.1396(1) 


1.637(37) 1.610(37) 
1.524(53) 1.518(61) 
1.454(87) 1.453(87) 


0.0325 


8 
16 


5.1456(6) 
5.1458(2) 


1.623(51) 1.607(49) 
1.535(69) 1.537(82) 


0.0350 


8 

12 

16 


5.1524(5) 
5.1508(5) 
5.1499(1) 


1.640(45) 1.625(44) 
1.664(55) 1.657(53) 
1.72(10) 1.72(10) 


0.0400 


12 

16 


5.1598(4) 
5.1593(5) 


1.896(63) 1.889(63) 
2.02(12) 2.01(12) 



Table 2: Critical couplings and fourth order cumulants 

Here we have determined the critical coupling /3c (^) from a linear interpolation of 
the pseudo-critical couplings on the 16^ x 4 lattice which are given in Tab. ^j. 

In addition to the Z(2) value for the Binder cumulant we show in Fig. ^ also 
the result for 3-dimensional 0(2) symmetric spin models. As the Binder cumulant 
depends quite sensitively on the underlying symmetry at the critical point the result 
found by us for the crossing point strongly suggests that the chiral critical point 
indeed belongs to the universality class of the 3-dimensional Ising model. The ratio 
of critical exponents 'j/u on the other hand is not sensitive to the universality class. 
As expected one finds, however, that all three susceptibilities show identical finite 
size scaling behaviour; for the quark mass closest to the critical point, m = 0.035, 
the ratio of susceptibilities calculated on lattices of size N^^ = 12 and 16 takes 
on the value 1.69(23), 1.66(22) and 1.66(24) for chiral, Polyakov-loop and action 
susceptibilities, respectively. This corresponds to a ratio of critical exponents 7/// = 
1.8(5) which is consistent with the 3-dimensional Ising as well as 0(2) and 0(4) 
values (7/2^ — 1.96). 

As discussed above the intersection point of Binder cumulants constructed from 
M{x) according to Eq. |^ will be independent of the choice of x in the infinite 
volume limit. In fact, from Eq. ^ it follows that for x = s the singular part of the 
free energy will lead to a unique crossing point, independent of the lattice size, i.e. 
finite volume effects are minimized if we manage to select the correct ordering-field 
like operator for constructing the Binder cumulant. For all other choices a volume 
dependence of the intersection points calculated from Binder cumulants on different 
size lattices will result from the scaling behaviour of /s(t, ^^nd they will converge 
to the universal value only in the infinite volume limit. It follows from Eq. ^jthat the 
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Figure 2: The Binder 



cumulant of the chiral condensate, -84(0) 



volume dependence is quadratic in A = x — s, i.e. in deviations from the optimal 
choice of x. This behaviour is confirmed through the analysis of the x-dependence of 
the intersection point of the Binder cumulant calculated on lattices of two different 
sizes. In Fig. ^ we show the intersection points of Binder cumulants calculated on 
lattices of size 8^ x 4 and 16^ x 4. As can be seen the volume dependence is small 
for a wide range of x values and the extremum, which corresponds to the optimal 
choice a; = s, is closest to the 3-d Ising value. In fact, determining the x-value 
where the extremum is reached provides a finite volume estimate for the mixing 
parameter s. From a jackknife analysis one finds for the location of the minimum 
Sinin = 0.430(23). 



4 The mixing parameters r and s 



Having located the chiral critical point we now can determine the mixing parameters 
r and s. This will allow us to construct the energy- and ordering-field like operators 
and obtain further evidence that the chiral critical point belongs to the universality 
class of the 3- dimensional Ising model. 

The parameters r and s are fixed by demanding that Ai should obey basic 
properties of the order parameter for spontaneous symmetry breaking at the critical 
point. In the symmetric phase the order parameter should stay constant along the 
line of vanishing external field {^ = 0). Of course, this should only hold in the 
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-0.8-0.6-0.4-0.2 0.2 0.4 0.6 0.8 



Figure 3: The intersection point of Binder cumulants B^lx) calculated on lattices 
of size N^ X 4 with Ncr = 8 and 16. The horizontal line shows the universal Ising 
value. 



vicinity of the critical point. We thus demand 

= , (10) 



m. 



:0+,5 = 



or equivalently 

{6S6M) = . (11) 

Using these relations together with Eqs. |] and |^ we obtain two conditions for the 
mixing parameters r and s in terms of the parameter B and expectation values of 

Sq and ijjip, 



-B , 

{djjjjdSc) - B{{5i)i)f) 
{{5SGf) - B{5^Pij 6Sg) 



'12) 



The parameter B controls the mixing of the gauge coupling and bare quark 
masses needed to define lines of constant ^. In the symmetry broken phase one line 
of constant ^ is known to us; the line of first order phase transitions defines the 
zero external field line (^ = 0) of the effective Hamiltonian. We thus can extract B 
from the quark mass dependence of the pseudo-critical couplings. Using Eq. ^ one 
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obtains 



i'-'^-l^l, . (13) 



Knowing B we also know r and can construct s using Eq. [T^. 

In order to determine the mixing parameter r from Eq. [1^ we have to approximate 
the derivative by finite differences as we can perform calculations of l3c{m) only at 
a discrete set of quark mass values. A first estimate may be given using our data 
on the largest lattice (16^ x 4). The slope of (3pc{m) defines r~^. We estimate this 
from a straight line fit to the values given in Tab. ^ which gives, r = 0.51(2). The 
mixing parameter r is large and definitely non-zero; as expected, a mixture of two 
operators is needed to construct an energy like observable in which the otherwise 
dominant ordering-field like contributions to Sq and ipip cancel. The importance 
of choosing the correct mixing parameter r becomes apparent from an analysis of 
joint probability distributions for 5£ and 5M.. These are little affected by changes 
of s but they strongly depend on r. In Fig. |^ we show contour plots for the joint 
probability distributions and various values of r and s. For r = s = these are just 
the fluctuations in 5Sg and dipip which are strongly correlated. Only for larger values 
of r fluctuations in the order parameter are independent of changes in the energy- 
like observable. In fact, we have used this criterion to improve our determination 
of r over our previous estimate based on the slope of Pdrh) at rh. We demand 
that fluctuations in 6Ai vanish for any fixed value of 6S. This generalizes Eq. ^ 
and maximizes the Z{2) symmetry of the contour plots shown in Fig.^ around the 
6Ai = axis. As an estimate for the mixing parameter r one obtains in this 
way r = 0.550(7). A determination of the parameter s from Eq. [1^ then yields 
s = 0.41(51). Although this value has large errors it is consistent with the result 
found from the extremum of intersection point of the Binder cumulant shown in 
Fig. ^. As a best estimate of the mixing parameters we therefore obtain 

r = 0.550 ±0.007 , s = 0.430 ± 0.023 . (14) 



5 Physical scale at the chiral critical point 

The analysis presented in the previous sections confirms that for finite values of the 
three light quark masses the chiral critical line in QCD belongs to the universality 
class of the 3-dimensional Ising model^. For the case of three degenerate quark 

'^This line may end in a tricritical point for some finite value of the strange quark mass ||5|. 
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Figure 4: Joint probability distributions of the fluctuations in the ordering-fleld hke 
and energy hke operators, which have been contructed according to Eq. |^. 



12 



m 


iV. 


Ppc 


^4(0) 


F-i((#^)2) 


0.005 


12 


3.2443(8) 


1.31(12) 


54.8(5.6) 


0.01 


16 


3.2778(8) 


2.14(10) 


11.5(1.0) 



Table 3: Critical couplings, Binder cumulant of the chiral condensate calculated and 
cliiral susceptibility calculated on N^ x 4 lattices with the p4-action [l^, |20l . 



masses we have located the critical point as given in Eq. ^. In order to determine 
a physical scale for this endpoint we have calculated hadron masses on a 16^ lattice 



at {/3c 



m] 



For the pseudo-scalar and vector meson mass we find m^s = 0.463(1) 



and m^ = 1.387(38), respectively. Expressing the pseudo-scalar mass in units of 
the critical temperature, rrips/Tc = 1.853(1), and using estimates for the critical 
temperature in 2 and 3-fiavour QCD [|T^ we thus estimate for the pseudo-scalar 
meson mass at the chiral critical point rrips — 290 MeV. 

The entire analysis of the chiral critical point discussed so far has been performed 
with unimproved gauge and fermion actions on rather coarse lattices. Improved 
actions are not expected to modify the results on the universal properties of the 
chiral critical point, which have been presented above. They may, however, well 
infiuence the quantitative determination of the chiral critical point. It has been 
found in studies of the first order deconfinement transition occurring in the pure 
gauge sector that the gap in physical observables like the latent heat or surface 
tension is cut-off dependent. Improved actions generally lead to smaller gaps and a 
reduced cut-off dependence of these observables |l^ . One thus may expect that also 
in the region of first order chiral transitions the gap in the chiral condensate gets 
reduced when calculated with improved actions. This will shift the critical point to 
smaller values of the pseudo-scalar meson mass. We therefore have investigated the 
chiral critical point also in calculations with improved gauge and staggered fermion 
actions (p4-action), which we have used previously for studies of the thermodynamics 
of two and three fiavour QCD [|l^, |l^]. In these calculations no evidence for a first 
order transition has been found down to bare quark masses m = 0.01 [0. We 
now have extended these calculations to smaller quark masses and made use of 
the universal properties of the Binder cumulants discussed above. One can locate 
the chiral critical point quite accurately through a calculation of Binder cumulants 
on finite lattices and there is no need to accurately determine the correct order 
parameter for such an analysis; i.e. the analysis can be performed with Binder 
cumulants constructed from the chiral condensate. We have performed calculations 
on a lattice of size 12^ x 4 with a bare quark mass m = 0.005 and on a lattice of 
size 16'^ X 4 with m = 0.01. Again a Ferrenberg-Swendsen reweighting is used to 
determine the pseudo-critical couplings and Binder cumulants at these couplings. 
The results are summarized in Tab. Bl 
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From the results given in Tab. ^ it is obvious that the smaller quark mass leads 
to a Binder cumulant below the value for the 3-dimensional Ising model and thus is 
in the first order region of the 3-flavour phase diagram. This also is confirmed by 
the large value of the chiral susceptibility. For the larger quark mass, on the other 
hand, the susceptibility stays small and the Binder cumulant is significantly above 
the Ising value. This quark mass still lies in the crossover region. We thus have 
obtained an upper and lower limit for the chiral critical point, suggesting a critical 
bare quark mass of m = 0.0075(25). Extrapolating the meson masses calculated 
in |]l^ to this value of the bare quark mass, we estimate for pseudo-scalar meson 



mass at the chiral critical point nips — 192(25) MeV. The critical mass thus is 
considerably smaller than estimated from calculations with unimproved gauge and 
fermion actions. 



6 Conclusions 

Through an analysis of Binder cumulants we have verified that the chiral critical 
point in three flavour QCD belongs to the universality class of the three dimensional 
Ising model. The analysis of joint probability distributions provides a powerful tool 
to construct the order parameter at the chiral critical point as well as the energy like 
scaling field. Although the chiral condensate itself is not the order parameter at this 
critical point, we explicitly have verified that Binder cumulants constructed from it 
are little influenced by finite volume effects and are good observables to locate the 
critical point as well as the universality class. 

Having determined the universality class of the chiral critical point for three 
degenerate quark masses one could use this information to determine the critical 
parameters also for non-degenerate quarks from calculations of Binder cumulants 
on finite lattices. The quality of the straight line fits shown in Fig. ^ suggests that a 
first order Taylor expansion of -64(0) in terms of degenerate up/down quark masses 
rriu^d and a strange quark mass rris around the three flavour critical point m = 0.033 
might be possible. As the chiral critical line corresponds to those sets of quark 
masses where the Binder cumulant attains the 3-d Ising value one has to determine 
the line on which -84(0) stays constant in the {mu,d, ms)-plane. To leading order this 
line is given by 

ms=rh-2 {111^4 - m) . (15) 

For instance, for rUu^ = 0.025 this estimate suggests that the region of first order 
chiral transitions ends at a critical value of the strange quark mass rris ~ 0.049. This 
is consistent with the result of 0] where no sign for a first order transition has been 
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found for {rriu^^ms) = (0.025,0.1) and not in contradiction with |jTl| who reported 
on a "weak first order hke behaviour" at (m^^^, m^) = (0.025,0.05) as opposed to 
two state signals at lower tUs values. 

The chiral critical point has been determined in calculations with unimproved 
as well as improved gauge and staggered fermion actions on lattices with temporal 
extent N^ = 4. The physical scale extracted from calculations of the pseudo-scalar 
meson mass at these endpoints is quite different in both cases. This indicates that 
cut-off effects are still significant and calculations closer to the continuum limit are 
definitely needed to fix a physical scale for the location of the chiral critical point. 
The present indication is, however, that improved actions lead to smaller values for 
the critical pseudo-scalar meson mass. This makes it increasingly unlikely that the 
transition in the physically realized case of two light and a heavier strange quark 
lies in the first order region of the QCD phase diagram. 
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